Superlinear and Quadratic Convergence of Riemannian Interior Point Methods

We extend the classical primal-dual interior point algorithms from the Euclidean setting to the Riemannian one. Our method, named the Riemannian interior point (RIP) method, is for solving Riemannian constrained optimization problems. Under the standard assumptions in the Riemannian setting, we establish locally superlinear, quadratic convergence for the Newton version of RIP and locally linear, superlinear convergence for the quasi-Newton version. These are generalizations of the classical local convergence theory of primal-dual interior point algorithms for nonlinear programming proposed by El-Bakry et al. and Yamashita et al. in 1996.Comment: 34page

A New Extension of Chubanov's Method to Symmetric Cones

We propose a new variant of Chubanov's method for solving the feasibility problem over the symmetric cone by extending Roos's method (2018) of solving the feasibility problem over the nonnegative orthant. The proposed method considers a feasibility problem associated with a norm induced by the maximum eigenvalue of an element and uses a rescaling focusing on the upper bound for the sum of eigenvalues of any feasible solution to the problem. Its computational bound is (i) equivalent to that of Roos's original method (2018) and superior to that of Louren\c{c}o et al.'s method (2019) when the symmetric cone is the nonnegative orthant, (ii) superior to that of Louren\c{c}o et al.'s method (2019) when the symmetric cone is a Cartesian product of second-order cones, (iii) equivalent to that of Louren\c{c}o et al.'s method (2019) when the symmetric cone is the simple positive semidefinite cone, and (iv) superior to that of Pena and Soheili's method (2017) for any simple symmetric cones under the feasibility assumption of the problem imposed in Pena and Soheili's method (2017). We also conduct numerical experiments that compare the performance of our method with existing methods by generating instances in three types: strongly (but ill-conditioned) feasible instances, weakly feasible instances, and infeasible instances. For any of these instances, the proposed method is rather more efficient than the existing methods in terms of accuracy and execution time.Comment: 44 pages; Department of Policy and Planning Sciences Discussion Paper Series No. 1378, University of Tsukub

A complexity bound of a predictor-corrector smoothing method using CHKS-functions for monotone LCP

We propose a new smoothing method using CHKS-functions for solving linear complementarity problems. While the algorithm in [6] uses a quite large neighborhood, our algorithm generates a sequence in a relatively narrow neighborhood and employs predictor and corrector steps at each iteration. A complexity bound for the method is also provided under the assumption that the problem is monotone and has a feasibleinterior point. As a result, the bound can be improved compared to the one in [6].Includes bibliographical references (p. 16-17

LP-based Tractable Subcones of the Semidefinite Plus Nonnegative Cone

Tractable Subcones and LP-based Algorithms for Testing Copositivit